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Abstract 

We evaluate the Hadamard function and the vacuum expectation value (VEV) of the current 
density for a charged scalar field, induced by flat boundaries in spacetimes with an arbitrary number 
of toroidally compactified spatial dimensions. The field operator obeys the Robin conditions on 
the boundaries and quasiperiodicity conditions with general phases along compact dimensions. In 
addition, the presence of a constant gauge field is assumed. The latter induces Aharonov-Bohm-type 
effect on the VEVs. There is a region in the space of the parameters in Robin boundary conditions 
where the vacuum state becomes unstable. The stability condition depends on the lengths of 
compact dimensions and is less restrictive than that for background with trivial topology. The 
vacuum current density is a periodic function of the magnetic flux, enclosed by compact dimensions, 
with the period equal to the flux quantum. It is explicitly decomposed into the boundary-free and 
boundary-induced contributions. In sharp contrast to the VEVs of the field squared and the 
energy-momentum tensor, the current density does not contain surface divergences. Moreover, for 
Dirichlet condition it vanishes on the boundaries. The normal derivative of the current density 
on the boundaries vanish for both Dirichlet and Neumann conditions and is nonzero for general 
Robin conditions. When the separation between the plates is smaller than other length scales, the 
behavior of the current density is essentially different for non-Neumann and Neumann boundary 
conditions. In the former case, the total current density in the region between the plates tends to 
zero. For Neumann boundary condition on both plates, the current density is dominated by the 
interference part and is inversely proportional to the separation. 


PACS numbers: 03.70.+k, ll.10.Kk, 04.20.Gz 


1 Introduction 

In a number of physical problems one needs to consider the model in the background of manifolds with 
boundaries on which the dynamical variables obey some prescribed boundary conditions. In quantum 
field theory, the imposition of boundary conditions on the field operator gives rise to a number of 
physical consequences. The Casimir effect is among the most interesting phenomena of this kind (for 
reviews see CD- It arises due to the modification of the quantum fluctuations of a field by boundary 
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conditions and plays an important role in different fields of physics, from microworld to cosmology. 
The boundary conditions in the Casimir effect may have different physical natures and can be divided 
into two main classes. In the first one, the constraints are induced by the presence of boundaries, 
like macroscopic bodies in QED, interfaces separating different phases of a physical system, extended 
topological defects, horizons in gravitational physics, branes in high-energy theories with extra di¬ 
mensions and in string theories. In the corresponding models the field operator obeys the boundary 
condition on some spacelike surfaces (static or dynamical). The original problem with two conducting 
plates, discussed by Casimir in 1948 [2], belongs to this class. Since the original research by Casimir, 
many theoretical and experimental works have been done on this problem for various types of bulk and 
boundary geometries. Different methods have been developed including direct mode-summation and 
the zeta function techniques, semiclassical methods, the optical approach, worldline numerics, the path 
integral approach, methods based on scattering theory, and numerical methods based on evaluation of 
the stress tensor via the fluctuation-dissipation theorem. The recent high precision measurements of 
the Casimir force allow for an accurate comparison between the experimental results and theoretical 
predictions. 

In the second class, the boundary conditions on the field operator are induced by the nontrivial 
topology of the space. The changes in the properties of the vacuum state generated by this type of 
conditions are referred to as the topological Casimir effect. The importance of this effect is motivated 
by that the presence of compact dimensions is an inherent feature in many high-energy theories of 
fundamental physics, in cosmology and in condensed matter physics. In particular, supergravity and 
superstring theories are formulated in spacetimes having extra compact dimensions. The compactified 
higher-dimensional models provide a possibility for the unification of known interactions. Models of 
a compact universe with nontrivial topology may also play an important role by providing proper 
initial conditions for inflation in the early stages of the Universe expansion [3|. In condensed matter 
physics, a number of planar systems in the low-energy sector are described by an effective field theory. 
The compactification of these systems leads to the change in the ground state energy which is the 
analog of the topological Casimir effect. A well-known example of this type of systems is a graphene 
sheet. In the long wavelength limit, the dynamics of the quasiparticles for the electronic subsystem is 
described in terms of the Dirac-like theory in two-dimensional space (see Ref. m- The corresponding 
effective 3-dimensional relativistic field theory, in addition to Dirac fermions, involves scalar and gauge 
fields (see [5] and references therein). The single-walled carbon nanotubes are generated by rolling 
up a graphene sheet to form a cylinder and for the corresponding Dirac model one has the spatial 
topology R 1 x S' 1 . For another class of graphene-made structures, called toroidal carbon nanotubes, 
the background topology is a 2-dimensional torus, T 2 . 

Many authors have investigated the Casimir energies and stresses associated with the presence 
of compact dimensions (for reviews see Refs. [U [6, Uj). In higher-dimensional models the Casimir 
energy of bulk fields induces an effective potential for the compactification radius. This has been used 
as a stabilization mechanism for the corresponding moduli fields and as a source for dynamical com¬ 
pactification of the extra dimensions during the cosmological evolution. The Casimir effect has also 
been considered as a possible origin for the dark energy in both Kaluza-Klein-type and braneworld 
models [8|. Extra-dimensional theories with low-energy compactification scale predict Yukawa-type 
corrections to Newton’s gravitational law and the measurements of the Casimir forces between macro¬ 
scopic bodies provide a sensitive test for constraining the parameters of the corresponding long-range 
interactions [9]. The influence of extra compactified dimensions on the Casimir effect in the classical 
configuration of two parallel plates has been recently discussed for scalar CD!, electromagnetic m 
and fermionic [12] fields. 

The vast majority of the works on the influence of the copmactification on the properties of the 
quantum vacuum in the Casimir effect has been concerned with global quantities such as the force 
or the total energy. More detailed information on the vacuum fluctuations is contained in the local 
characteristics. Among the most important local quantities, because of their close connection with the 
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structure of spacetime, are the vacuum expectation values (VEVs) of the vacuum energy density and 
stresses. For charged fields, another important characteristic is the VEV of the current density. Due to 
the global nature of the vacuum, this VEV carries information on both global and local properties of 
the vacuum state. Besides, the VEV of the current density appears as a source of the electromagnetic 
field in semiclassical Maxwell equations, and, hence, it is needed in modeling a self-consistent dynamics 
involving the electromagnetic field. 

In models with nontrivial topology, the nonzero current densities in the vacuum state may appear as 
a consequence of quasiperiodicity conditions along compact dimensions or by the presence of gauge field 
fluxes enclosed by these dimensions. Note that the gauge field fluxes in higher-dimensional models will 
also generate a potential for moduli fields and this provides another mechanism for moduli stabilization 
(for a review see [13]). The VEV of the fermionic current density in spaces with toroidally compactified 
dimensions has been considered in p3]. In the special case of a 2-dimensional space, application are 
given to the electrons in cylindrical and toroidal carbon nanotubes, described within the framework 
of the effective field theory in terms of Dirac fermions. The vacuum currents for charged fields in 
de Sitter and anti-de Sitter spacetimes with toroidally compact spatial dimensions are investigated in 
mills]. Finite temperature effects on the charge density and on the current densities along compact 
dimensions have been discussed in m and [TEJ f° r scalar and fermionic fields, respectively. The 
changes in the fermionic vacuum currents induced by the presence of parallel plane boundaries, with 
the bag boundary conditions on them, are investigated in 191- 

In the present paper we consider the effect of two parallel plane boundaries on the vacuum ex¬ 
pectation value of the current density for a charged scalar field in background spacetime with spatial 
topology R p+l x T q , where T q stands for a q-dimensional torus. The organization of the paper is as 
follows. In the next section the geometry of the problem is described and the Hadamard function 
is evaluated in the region between the plates for general Robin boundary conditions. By using the 
expression for the Hadamard function, in Section [3] we evaluate the current density in the geometry 
of a single plate. The corresponding asymptotics are discussed in various limiting cases and numerical 
results are presented. In Section |4] the current density is investigated in the region between two plates. 
The main results of the paper are summarized in Section [5] An alternative representation of the 
Hadamard function is given in Appendix. 

2 Formulation of the problem and the Hadamard function 

We consider (D + l)-dimensional flat spacetime with spatial topology R p+1 x T q , p + q + 1 = D (for 
a review of quantum field-theoretical effects in toroidal topology see Ref. ED- The set of Cartesian 
coordinates in the subspace R p+1 will be denoted by x p+ i = (x 1 ,..., x p+1 ) and the corresponding 
coordinates on the torus by x g = ( x p+2 , ...,x D ). If L[ is the length of the Zth compact dimension then 
one has — oo < x l < oo for l = 1, ..,p, and 0 ^ x l ^ Li for l = p + 2,..., D. Our main interest in this 
paper is the VEV of the current density for a quantum scalar field <p(x) with the mass m and charge 
e. The equation for the field operator reads 

(g^D^D u + m 2 )ip = 0, (2.1) 

where g pu = diag(l, —1,... , —1), D M + ieA M and A M is the vector potential for a classical gauge 

field. We assume the presence of two parallel flat boundaries placed at x p+l = oq and x p+1 = 02 , on 
which the field obeys Robin boundary conditions 

(1 + (3jrij D^tp^x) = 0, x p+1 = z = cij , (2.2) 

with constant coefficients (3j, j = 1, 2, and with rij being the inward pointing normal to the boundary 
at x p+l = cij. Here, for the further convenience we have introduced a special notation z = x p+l for the 

1 In analogy with the standard Casimir effect, in the discussion below we will refer the boundaries as plates. 
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(p+ l)th spatial dimension. Note that Robin boundary conditions in the form (12.21) are gauge invariant 
(for the discussion of various types of gauge invariant boundary conditions see [20]). In what follows 
we will consider the region between the plates, a\ ^ z ^ 02 - For this region one has rij = (— l)- 7_1 <5p +1 . 
The expressions for the VEVs in the regions z ^ a\ and z ^ ci 2 are obtained by the limiting transitions. 
The results for Dirichlet and Neumann boundary conditions are obtained from those for the condition 
(12.21) in the limits f3j —>• 0 and (3j —>• 00 , A M = 0, respectively. Robin type conditions appear in a 
variety of situations, including the considerations of vacuum effects for a confined charged scalar field 
in external fields [21], gauge field theories, quantum gravity and supergravity [201122], braneworld 
models [23] and in a class of models with boundaries separating the spatial regions with different 
gravitational backgrounds |Mj- In some geometries, these conditions may be useful for depicting 
the finite penetration of the field into the boundary with the ” skin-depth” parameter related to the 
coefficient j3j. It is interesting to note that the quantum scalar field constrained by Robin condition 
on the boundary of cavity violates the Bekenstein’s entropy-to-energy bound near certain points in 
the space of the parameter j3j [25]. 

In addition to the boundary conditions on the plates, for the theory to be completely defined, 
we should also specify the periodicity conditions along the compact dimensions. Different conditions 
correspond to topologically inequivalent field configurations [26]. Here, we consider generic quasiperi¬ 
odicity conditions, 


<f{t, x 1 ,..., x l + L h ..., x D ) = e iai v(t, x 1 ,..., x l ,..., x D ), (2.3) 

with constant phases cq, l = p + 2,..., D. The special cases of the condition (|2.3I) with cq = 0 
and on = it correspond to the most frequently discussed cases of untwisted and twisted scalar fields, 
respectively. As it will be seen below, one of the effects of nontrivial phases in (12.3jl is the appearance 
of nonzero vacuum currents along compact dimensions (for a discussion of physical effects of phases 
in periodicity conditions along compact dimensions see m and references therein). 

For a scalar field, the operator of the current density is given by the expression 

= ie[ip + (x)D^(x) - (D^tp(x)) + ip(x)}, (2.4) 

Z = 0,1,...,.D. Its VEV is obtained from the Hadamard function 

G(x,x') = (0\if(x)ip + (x r ) + 93 + (x , )<^(x)|0), (2.5) 

with |0) being the vacuum state, by using the formula 

(0|j M (x)|0) = 0'/*(®)) = I™ (<9 m - d 1 + 2ieA^)G(x,x'). (2.6) 

In the discussion below we will assume a constant gauge held A^. Though the corresponding 
held strength vanishes, the nontrivial topology of the background spacetime leads to the Aharonov- 
Bohm-like effects on physical observables. In the case of a constant gauge held A^, the latter can 
be excluded from the held equation and from the expression for the VEV of the current density by 
the gauge transformation A^ = A'^ + d^Xi v{ x ) = e~ iex ip'{x), with the function x = A^x^ 1 . In the 
new gauge one has A' = 0. However, unlike to the case of trivial topology, here the constant vector 
potential does not completely disappear from the problem. It appears in the periodicity conditions 
for the new held operator: 

(/(t,® 1 , ...,x l + L h .. .,x D ) = e lQ V(M\ -..,x 1 ,.. -,x D ), (2.7) 

where now the phases are given by the expression 

on = on + eAiL t . (2.8) 
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In the discussion below we shall consider the problem in the gauge (ip'(x), A' = 0) omitting the prime. 
For this gauge, in (12.11) . (12.21) . (12.41) one has D p = d p and in the expressions (12.61) the term with the 
vector potential is absent. 

From the discussion above it follows that in the problem at hand the presence of a constant gauge 
field is equivalent to the shift in the phases of the periodicity conditions along compact dimensions. 
The shift in the phase is expressed in terms of the magnetic flux <1^ enclosed by the Ith. compact 
dimension as 

eAiLi = —eAiLi = -2 tt^i/^ 0 , (2.9) 

where <ho = 27r/e is the flux quantum and A; is the Zth component of the spatial vector A = 
(— Ai,..., — Ad)- In the discussion below the physical effects of a constant gauge field will appear 
through the phases oq. In particular, the VEVs of physical observables are periodic functions of these 
phases with the period 27r. In terms of the magnetic flux, this corresponds to the periodicity of the 
VEVs, as functions of the magnetic flux, with the period equal to the flux quantum. 

For the evaluation of the Hadamard function in (|2.6D we shall use the mode-sum formula 

G(x,x') (2.10) 

k -s=± 

where (x) form a complete set of normalised positive- and negative-energy solutions to the classical 
field equation obeying the boundary conditions of the model. In the region between the plates, 
introducing the wave vectors k p = (Aq,. .., k p ) and k g = (k p+ 2 , ■ • •, ko), these mode functions can be 
written in the form 


V\t\ x ) = Ck cos [fcp +1 (z - aj) + Jj(k p+ i)] e lk ir x ll=F^kt ) 


( 2 . 11 ) 


where ky = (k p ,k g ), k = (kp, k p+ \, k 9 ), = Vk 2 + to 2 , and xy stands for the coordinates parallel 

to the plates. For the momentum components along the dimensions x l , i = 1,... ,p, one has —00 < 
ki < + 00 , whereas the components along the compact dimensions are quantized by the periodicity 
conditions m-- 

ki = (2irni + ai)/Li, ra/= 0, ±1, ±2,...(2.12) 

with l = p + 2,..., D. We will denote by cuo the smallest value for the energy in the compact subspace, 
y^k^ + m 2 ^ wo- Assuming that |cq| ^ 7r, we have 

= \jYl D l=p+ 2^ /L2 l +m2 - (2 - 13) 

This quantity can be considered as the effective mass for the field quanta. 

Now we should impose on the modes (12.111) the boundary conditions (|2.2I) with D p = d p . From 
the boundary condition on the plate at z = a,j, for the function 7j(^' p +i) in (12.111) one gets 


= 2*7,-(fe p+ i) _ ik P +iPj( 1) J + 1 

ik p+ i/3j(-iy - 1' 


(2.14) 


From the boundary condition on the second plate it follows that the eigenvalues for k p+ 1 are solutions 
of the equation 

1 - 1 - i /in7/I -U ih-t 11 

(2.15) 


e 2 i. y = 1 + ib 2 y 1 + ibiy 
1 - ib2y 1 — ibiy' 


where 


V = k p+1 a, bj = f3j/a , 


(2.16) 
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and a = 02 — aq is the separation between the plates. Formula ([2.151) can also be written in the form 

(l - b\b 2 y 2 ) sin y - (b 2 + 6 i)y cosy = 0. (2.17) 


Unlike to the cases of Dirichlet and Neumann conditions, for Robin boundary condition the eigenvalues 
of k p .|_i are given implicitly, as solutions of the transcendental equation (12.171) . This equation has 
an infinite number of positive roots which will be denoted by y = X n , n = 1 , 2 ,..., and for the 
corresponding eigenvalues of k p+ 1 one has k p+ \ = A n /a. For bj ^ 0 or {b\ +b 2 ^ 1, b\b 2 ^ 0} there are 
no other roots in the right-half plane of a complex variable y, R ey ^ 0 (see [2BJ) - In the remaining 
region of the plane ( 61 , 62 )) the equation (12.171) has purely imaginary roots ±iyi, yi > 0. Depending 
on the values of bj, the number of yi can be one or two. In the presence of purely imaginary roots, 
under the condition wo < yi, there are modes of the field for which the energy Wk becomes imaginary. 
This would lead to the instability of the vacuum state. In the discussion below we will assume that 
wo > yi- Note that in the corresponding problem on background of spacetime with trivial topology 
the stability condition is written as m > y\. Now, by taking into account that ujq > m, we conclude 
that the compactification, in general, enlarges the stability range in the space of parameters of Robin 
boundary conditions. 

The coefficient Ck in (|2.11|) is found from the orthonormalization condition 

J d D X<p£\x)<p§ ] *(x) = ^-<K k P - K^ 6 nn' S n p+2 ,n' p+2 -A D ,n' D , ( 2 -!8) 


where the integration over x p+1 goes in the region between the plates. 
( 12 . 111 ) . one gets 


, r ,2 = {l + cos[y + 27 j(j/)]sin(j/)/y} 1 
1 k| (2vr)PaU g w k 


Substituting the functions 


(2.19) 


where y is a root of the equation (|2.171) and V q = L p+ i....Ld is the volume of the compact subspace. 
The function Jj(y) is defined by the relation 


e 2 *R (y) = Hah - 1 . ( 2 . 20 ) 

iybj + 1 

First we shall consider the case when all the roots of (|2.17l) are real and y = \ n . 

Having the complete set of normalized mode functions, the mode-sum (12.101) for the Hadamard 
function is written in the form 


G(x, x') 



X n cos(wkAf)e ikp ' Axp+ * k9 ' Ax,J 
An T COS [An T 2 7 j(An)] sin An 


( 2 . 21 ) 


where Ax p = x p —x p , Ax g = x q — x' q , At = t — t', and n g = (n p+2 ,... —00 < ni < + 00 . In (12.211) . 

the energy for the mode with a given k is written as 


w k = + XI/cl 2 + 


and 


L0 r 


D 


= y/k2 + m 2 , k 2 q = J2 ( 

l=p+2 


( 2iTill + Oil 


L, 


Here and in what follows we use the notation 


9 j(z,z',u) = cos (uAz) + - ^ 


I r siy\z+z’-2a i \ m Pj ~ g 


s==tl 


iu/3j + S 


( 2 . 22 ) 


(2.23) 


(2.24) 
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Note that gj(z,z',—y) = gj(z,z',y) and gj(z,z', 0) = 0. 

In (12.211) . the eigenvalues X n are given implicitly and this expression is not convenient for the 
evaluation of the VEVs. In order to obtain an expression in which the explicit knowledge of X n is not 
required, we apply to the series over n the Abel-Plana-type summation formula |28l 29] 


V vrA n /(A n ) = tt/(0)/2 

“ Xn + cos[A n + 27 j (A n )] sin A n 1 -b -2 -b] ,/o 

r du f(iu) - f(-iu) 
Jo Ci(u)c 2 (u)e^ ~ 1 

where, for the further convenience, the notation 


(2.25) 


, . bjU — 1 . 

= 6m (2 - 26) 

is introduced. In (12.251) we have assumed that bj ^ 0. The changes in the evaluation procedure in the 
case bj > 0 will be discussed below. For the series in (I2.2ip . we take in the summation formula 

/(An) = C08{UkAt) g j (z,z , ,X n /a). (2.27) 

Note that /(0) = 0 and the first term in the right-hand side of ()2.25|) is absent. 

The use of the summation formula (|2.25j) with (|2.27D allows us to write the Hadamard function in 
the decomposed form 


G(x, x') 


Gj(x,x') + -^r [ / dugj(z,z',iu/a) 

KVqJ 1/7 T) n? J auj k|| 
gikp-Axp+ikg-Ax, cosh(A t u 2 / a 2 - £U k|| ) 


ci(u)c2(w)e 2u — 1 


u 2 - a 2 u k .. 


(2.28) 


where cu k|| 


k 2 + u 2 


Here, the part 


Gj(x, x') 


1 

^V q 



gikp ■ Ax p +ik 9 ■ Ax, 


L 


00 cos(w k At) , 

dk p+ 1 -£,'( 2,2 ,k p +i), 

w k 


(2.29) 


comes from the first integral in the right-hand side of (]2.25|) and corresponds to the Hadamard function 
in the geometry of a single plate at x p+l = aj when the second plate is absent. This function is further 
decomposed by taking into account that the part in (]2.24|) coming from the first term in the right-hand 
side of (12.241) . 


Gq(x, x') 


1 

y q 


dkp+i zk p+ i-Ax p+ i+zk g -Ax g cos(o; k At) 

(2it)p+ 1 ^ uj k 


(2.30) 


is the Hadamard function for the boundary-free geometry. After the integration over the components 
of the momentum along uncompactified dimensions, this function can be presented in the form 


G o(x,x') 


2V- 1 

( 2 vr)P / 2 + 1 


£ 


e *VAx, „ 


fp/2^n q y^lAXp+il 2 - (At) 2 ), 


(2.31) 
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with the notations 


fu(x) = K v (x)/x v , 


(2.32) 


where K v {x) is the Macdonald function. 

Consequently, the Hadamard function in the geometry of a single plate is written as 


Gj(x,x') = G 0 (x,x') + 


1 


cik r . 


2 irVq J (27 t)p 


E 


e ik p • Ax p +ik q • Ax g 


___ /*c 

x E / 


’ cos(w k At) _ sifcn+1 u+z /_ 2o ,.| *Vn& - s 
Uj k ik p+ it3j + s 


(2.33) 


where the second term in the right-hand side is induced by the presence of the plate at x p+1 = aj. For 
the further transformation of the boundary-induced part in (12.3311 we rotate the integration contour 
over fcp+i by the angle s7t/2. In the summation over s the integrals over the intervals (0, ±iwk||) cancel 
each other and we get 


Gj(x, x 1 ) 


G 0 {x,x') + — 


7rl4 



gikp ■ Axp+«k g ■ Ax, 




U Pj + 1 -u\z+z'-2ai\ 

u(3j - 1 


(2.34) 


This expression is well suited for the investigation of the current density. With the representation 
(|2.34p . the Hadamard function in the region between the plates, given by (|2.28p . is decomposed 
into the boundary-free, single plate-induced and second plate-induced contributions. An alternative 
expression for the Hadamard function is obtained in Appendix. 

In deriving (12.28(1 and (12.3411 we have assumed that (3j ^ 0. In the case f3j > 0, the quantum 
scalar held in the geometry of a single plate at z = aj has modes with fc p +i = i/fij for which the 
dependence on the coordinate x p+1 has the form In the case 1 /f3j > ujq, for a part of these 

modes the energy is imaginary and the vacuum is unstable. In order to have a stable vacuum, in 
what follows, for non-Dirichlet boundary conditions, we shall assume that 1 //3j < uq and the mode 
with k p . |_i = i//3j corresponds to a bound state. For /3j > 0 and in the absence of purely imaginary 
roots of (12.17p . in the right-hand side of the summation formula (12.2511 the residue terms at u = ±i/bj 
should be added (see [28]). Now the integrand in (|2.33p has a simple pole at k p +\ = is/f3j and after 
the rotation the contribution of the residue at that pole should be added. This contribution cancels 
the additional residue term in the right-hand side of ((2.2511 . In the case when the equation (12.171) has 
purely imaginary roots the corresponding contributions have to be added to the mode-sum (12.2111 for 
the Hadamard function. But the corresponding contributions should also be added in the left-hand 
side of (12.251) and the further evaluation procedure remains the same. Hence, the expressions (12.281) 
and ((2.3411 are valid for all values of the coefficients in the Robin boundary conditions. The only 
restrictions come from the stability of the vacuum state: l/(3j < loq and yi < ujq- In the presence of 
compact dimensions with cq ^ 0 one has ujq > m and these conditions are less restrictive than those 
in the case of trivial topology. 

The current density in the boundary-free geometry is obtained by using the Hadamard function 
(12.311) and has been investigated in [17j . The corresponding charge density and the current densities 
along uncompact dimensions vanish. As it can be seen from (12.2811 and (12.341) . the same holds in 
the case of the boundary-induced contributions in the VEVs. Hence, the only nonzero components 
correspond to the current density along compact dimensions. 
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3 Vacuum currents in the geometry of a single plate 


In this section we investigate the VEV of the vacuum current density in the geometry of a single plate 
at x p+1 = dj. This VEV is obtained with the help of the formula (12, 6 |i by using the Hadamard function 
from (12.341) . The component of the VEV of the current density along the Zth compact dimension is 
presented in the decomposed form 

(A = ii'h + (Z)?’. (3-i) 

where (j l )o is the current density in the boundary-free geometry and (j 1 )^ is the contribution induced 
by the presence of the plate. 

The current density in the boundary-free geometry has been investigated in m and for the com¬ 
pleteness we will recall the main results. The current density is given by the formula 


</> o 


4 eLim D+l ^ ~ \ 

( 27r) (D+1)/2 2^ sm {nm) 

' ni =1 

x ^ cos(n.q—i • a q -i)fD+i(mg nq (L q )), 

n q ~i 


(3.2) 


where « g _i = (a p+ 2 , ■ ■ ■ , d;_i, ai+i, ■ ■ ■, old), n g _i = (n p+ 2 , ■ ■ ■ ,ni-i,ni + i,... ,ud), and g nq (L q ) = 
( 52 iLp +2 The current density (j l )o is an odd periodic function of d/ with the period 27 t and 

an even periodic function of d r , r ^ l, with the same period. This corresponds to the periodicity in 
the magnetic flux with the period of flux quantum. An alternative expression for the current density 
in the boundary-free geometry is given by the formula 


/-i\ 4e Li/V q ^ sin (ndq) ^ , T \ 

~ (2vr)(P+ 3 )/ 2 ^ ( nLi)P+ 2 ^ W n,-i)> 

71 —1 Wq—\ 

(3.3) 

where we have defined the function 


9v{x) = x v K v {x), 

(3.4) 

and 


2 _ 2 r ,2 

U n q -l ~ W n, “ k l ■ 

(3.5) 


In the model with a single compact dimension (q = 1) the representations (|3.2D and (|3.3D are identical. 

When the length of the (th compact dimension, Li , is much larger than the other length scales, 
the behavior of the current density crucially depends whether the parameter 


uoi 



D 

i=p-\-2,^l 


a: 


'/Ll + nr/f 


1/2 


(3.6) 


is zero or not. For woz 
order we have 


0 , which is realised for a massless field with dj 

: i ~ 2eT((p + 3)/2) ^ sin (red;) 

°~ir( p+ 3 )/ 2 L p+ 1 V„ n p + 2 

L h n=l 


0 , i 7 ^ (, to the leading 


(3.7) 


In this case, the leading term in the expansion of V q (j l )o/Li coincides with the current density in 
(p + 2)-dimensional space with a single compact dimension of the length L[. For ujqi ^ 0 and for large 
values of Li one has 


if) o 


2eV q 1 sin(d;)u ;^ //2+1 
{2tt)p/ 2 + 1 L p 1 /2 


e ~Liui 0 i 


(3.8) 
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and the current density is exponentially suppressed. In the opposite limit of small values for L;, to 
the leading order we get 


, .n 2eY{{D + l)/2) sin(nch) 
{J) °~ vr^+ 1 )/ 2 L / 3 ^ 

( n=l 


(3.9) 


The leading term does not depend on the mass and on the lengths of the other compact dimensions 
and coincides with the current density for a massless scalar field in the space with topology R D ~ l x S' 1 . 

Now we turn to the investigation of the plate-induced contribution in the current density. By using 
the expression for the corresponding part in the Hadamard function from (12.341) . we get the following 
expression 


</>?> - 


cCp 

2 W„ 


___ r c 

E *>/ 

„ J OJr 


dy {y 2 - u>l 


)(p- 1 )/ 2 e ~ 2 y z : 
1 ' 


, yfij + 1 

vPj - 1 : 


(3.10) 


with the notations Zj = \z — Oj| for the distance from the plate and 


7r-(p+B/2 

= T((p + l)/2)‘ 


(3.11) 


Recall that, in order to have a stable vacuum state with (ip) = 0, we have assumed that l//3j < ujq. 
Under this condition, the integrand in (13.101) is regular everywhere in the integration range. The 
integral in (13.101) is evaluated in the special cases of Dirichlet and Neumann boundary conditions with 
the result 

ti l )j ) = T k l UJ n q fp/2(2 UJ nq Zj), (3-12) 

n <7 

where the upper and lower signs correspond to Dirichlet and Neumann boundary conditions, respec¬ 
tively. Note that, in the problem with a fermionic field, obeying the bag boundary condition on the 
plate, the boundary-induced contribution vanishes for a massless field )19j . 

Let us consider the behavior of the plate-induced contribution in asymptotic regions of the pa¬ 
rameters. At large distances from the plate, zj L l , one has Zjio nq 1. Assuming that |dj| < 7 r, 
the dominant contribution in (13.101) comes from the region near the lower limit of the integration and 
from the term with rii = 0, i = p + 2,..., D. To the leading order we find 

/.-Mi) eogJ^^e- 2 ^ upPj + 1 

° )j ~ (4f)W^ “ofr - 1 ’ 1 j 

and the current density is exponentially small. Note that the suppression is exponential for both 
massive and massless field. 

For points close to the plate, Zj <C in (13.101) the contribution of the terms with large values of |n,| 
dominates and this formula is not convenient for the asymptotic analysis and for numerical evaluations. 
In the case ^ 0, an alternative expression is obtained by using the representation (IA. 6 D for the 
Hadamard function. The first term in the right-hand side of this representation corresponds to the 
geometry with uncompactified Ith. dimension and does not contribute to the current density along that 
direction. In the geometry of a single plate at x p+1 = a,j the part in the Hadamard function induced 
by the compactihcation is given by the first term in the figure braces of (|A. 6 I) . From this part, by 
making use of (12.61) . for the VEV of the Ith. component of the current density we get 


z _ 2 1 p/2 eLi ^ sin (ncq) 

\J )j ~ _n/24-2TA / j 


nP ,2 + 2 Vq ^, {nLl y + i 


roo -- 

/ d y9(zj,y)g p/2+ i(nLiyJy 2 + uj (3.14) 
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where we have defined the function 


g(zj,y) = gj(z,z,y) = 1 + ^ ^ 

s=±l 


1 p 2siy Zj WPjJZJL 

iy/3j + s 


= 1 - 


(! - y 2 / 3 |) COS (2 yzj) + 2y/3j sin(2 yzj) 

1 + y 2 P] 


(3.15) 


The part with the first term in the right-side of (13.151) corresponds to the current density in the 
boundary-free geometry. In this part the integration over y is done with the help of the formula 


/ dy gz_ |_i (nL t y 2 + b 2 ) = ^/V/2{nLi) 1 q P +3 (nigh), 

Jo - 2 

and one gets the expression (13.31) . 

Extracing the boundary-free part, for the plate-induced contribution from (13.141) we find 


(3.16) 



2 P//2 eL; ^ sin {non) 
7rP /2 + 2 V q ^ [nL l )P+ l 



dygE + i(nLi 



e 2 siyzj 

s=± 1 


iyfy - s 

iy/3j + s' 


(3.17) 


In the case of single compact dimension one has q = 1, p = D — 2, and the corresponding formula for 
the plate-induced contribution in the current density is obtained from (I3.17P omitting the summation 
over n g _! and putting w nr) = m. 

An important issue in quantum field theory with boundaries is the appearance of surface diver¬ 
gences in the VEVs of local physical observables. Examples of the latter are the VEVs of the field 
squared and of the energy density. These divergences are a consequence of the oversimplification of 
a model where the physical interactions are replaced by the imposition of boundary conditions for 
all modes of a fluctuating quantum field. Of course, this is an idealization, as real physical systems 
cannot constrain all the modes (for a discussion of surface divergences and their physical interpre¬ 
tation see mm and references therein). The appearance of divergences in the VEVs of physical 
quantities indicates that a more realistic physical model should be employed for their evaluation on 
the boundaries. An important feature, which directly follows from the representation (13.171) . is that 
the VEV of the current density is finite on the plate. This is in sharp contrast with the behavior of the 
VEVs for the field squared and energy-momentum tensor. The finiteness of the current density on the 
boundary may be understood from general arguments. The divergences in local physical observables 
are determined by the local bulk and boundary geometries. If we consider the model with the topology 
RP +2 x T 9_1 with the Ith dimension having the topology R 1 , then in this model the Ith component of 
the current density vanishes by the symmetry. The compactification of the Ith dimension to S 1 does 
not change both the bulk end boundary local geometries and, hence, does not add new divergences to 
the VEVs compared with the model on R p+2 x T 9_1 . 

In deriving (|3.171) we have assumed that / 3j ^0. In the case [3j > 0 the contribution of the bound 
state should be added to (13.171) . For 1 /flj < uj 0 i, this contribution is obtained from the corresponding 
part in the Hadamard function, given by (IA.7I) . and has the form 


od? = - 


2 2 p / 2 eLie 2z i!di sin (nch) 


7T P/ 2+1 V q P 3 


sill y iLU-l J 


9t+i(»rV“».-i - 1 A ? )- I 3,18 ) 


In what follows for simplicity we shall consider the case /3 j ^0. Recall that, the representation (13.101) 
is valid for all values of (3j from the range of the vacuum stability. 

For Dirichlet and Neumann boundary conditions, after the evaluation of the integral in (13.171) by 
using the formula 


j r dy cos(2yZj)gR +1 (nLiky 2 + b 2 ) = y^( nLi) p+2 


Q p+3 (&y/4 z 2 + n 2 L 2 ) 
(4z? + n 2 L 2 )(P+ 3 )/ 2 ' 


(3.19) 
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one gets 


tf >? -C3.20) 


nsin (ndj) 


where the upper and lower signs correspond to Dirichlet and Neumann conditions, respectively. For 
a single compact dimension with the length L and with the phase d in the periodicity condition for a 
massless field this gives 


,_ 2T((D + l)/2)e nsin(nd) 

V'i -T ^{ 0 + 1 ) 12 L D zL ( n 2 + 4 2 / L 2 )(D+l )/2 ' 


(3.21) 


Now, combining the expressions (13.311 and (|3.20l) . we see that in the case of Dirichlet boundary condition 
the boundary-free and plate-induced parts of the current density cancel each other for Zj = 0 and, 
hence, the total current vanishes on the plate. For Neumann condition the current density on the 
plate is given by 



,z= aj 


2 (j l )o = 


8eLi/V q sin (ncq) 

(2tt)(p + 3 )/ 2 ( nL t )P+ 2 


Y grklinLiuJn^). 

n,_i 2 


(3.22) 


Note that the normal derivative of the current density on the plate vanishes for both Dirichlet and 
Neumann boundary conditions: (d z (j l )j) z = a] = 0. This is not the case for general Robin condition. 

Let us consider the behavior of the plate-induced contribution in the current density in the 
limit Li < L[. In this investigation it is more convenient to use the representation (|3.17l) . For 
J2iL p +2,^l the dominant contribution in the integral of (13.1711 comes from the region near the 

lower limit of the integration and from the term n = 1, rii = 0, i = p+ 2,..., D, in the summation. The 
argument of the function g p / 2 +i(x ) in the integrand is large and we can use the asymptotic expression 
g v (x) ~ y / 7 r/ 2 x i ' - 1 / 2 e -3: . After some intermediate calculations, for the leading term we get 

i) „ 2 e(l- 2 S op .) cug/ 2 + 1 sind; 

j (2ir)P/ 2+1 V q L^ 2 e D^oi(i+2zf/D 2 ) 

Here, we have additionally assumed that Li <C \(3j\ for /3j ^ 0. For dj = 0, i = p + 2,..., D, i / l, the 
dominant contribution in (13.1711 comes from the term n* = 0, i = p + 2,... ,D, with the leading term 




w /j^p+ixS 1 


4e 


E 


sin (nd;) 


f 

Jo 


xgE +1 (nLi^y 2 + m 2 ) Y ^ 


( 2 7 t)p/ 2+2 ^ (nL z )P+i J 0 

2 aivzH iyfc ~ S 
iy(Jj + s 


dy 


s=±l 


(3.24) 


Here, {j l )^ RP +i xS i is the plate-induced contribution in the current density for (p + 2)-dimensional 
space with topology R p+1 x S 1 (see (13.1711 for the case q = 1 and, hence, ix nq _ 1 = m). 

If the length of the ith compact dimension is large, i ^ l, the dominant contribution to the sum 
over m comes from large values of |nj| and in (13. 1 711 we can replace the summation over n, by the 
integration in accordance with 

00 t roo 

Y dx ( 3 - 25 ) 

ni=—oo ^ J 0 

The integral over x is evaluated by using the formula ()3.16jl . As a result, from (|3.17j) . to the leading 
order, we obtain the current density along the Ith compact dimension for the spatial topology R p+2 x 
T q ~ l with the lengths of the compact dimensions (L p+ 2 ,..., Lj_i, Lj+i,..., Lp). 


12 































Now let us consider the limiting case when Li is large compared with the other length scales in the 
problem, L; 3> Li,Zj, i ^ l. The dominant contribution in (13.171) comes from the term Uj = 0, i ^ l. 
For uqi 7 ^ 0 we find 


/.■Mi) ~ 2e ( 2 ^j,oo !) sm ay p/2+i Ll u ol 

U ~ (2vr)P/ 2 + 1 Vq Lf/ 2 01 


(3.26) 


where, for non-Neumann boundary conditions (/ 3j ^ oo), we have assumed that (3jUoi <C (Liloqi) 1 ^ 2 . 
For ojqi = 0 the leading term is given by the expression 



2 e ( 2 fy . |00 - 1 ) 
7 r (p+3)/2^Lf +1 


r ((p + 3)/2) J2 

n= 1 


sin (nag) 
nP + 2 


(3.27) 


Comparing with the corresponding asymptotics (13.71) and (13.81) . we see that for non-Neumann boundary 
conditions, in the both cases loqi / 0 and ujqi = 0 , the leading terms in the boundary-induced and 
boundary-free parts of the current density cancel each other. 

An equivalent representation for the plate-induced current density is obtained from (13.171) rotating 
the integration contour in the complex plane y by the angle 7r/2 for the term with s = 1 and by the 
angle —ir/2 for the term with s = —1. The integrals over the intervals (0, x ) are cancelled and 

we find 

/.•Mi) = 2~ v/2 eLi sin(ndz) f°° 

Xe ~ 2VZj V yp ■ - \ w p/ 2 +^ nLl \Jy 2 ~ (3.28) 

where 

w v (x) = x u J u (x), (3.29) 

and J u (x) is the Bessel function. The equivalence of the representations (13.101) and (13.28[) can also be 
directly seen by applying to the series over ni in (13.101) the relation 


+00 2-Zv °° /'OO 

y, hg{\ki\) =-y sin(naj) / dxxsm(nLix)g(x). 

ni=—o o n=1 ' ^ 


(3.30) 


The latter is a direct consequence of the Poisson’s resummation formula. After using (13.30[) in (13.10|) . 
we introduce a new integration variable u = ^ Jy 2 — x 2 — w 2 ? _ 1 and then pass to polar coordinates in 

the (?i, x)-plane. The integration over the polar angle is expressed in terms of the Bessel function and 
the representation (13.281) is obtained. 

Another expression is obtained by applying to the series over m in (13.101) the summation formula 
(1A.1I) . For the series in (13.101) one has g(u) = u and the first integral vanishes. As a result, the 
plate-induced part in the VEV of the current density is presented as 


{J l )f = 


eC p Li sin ay 


2PttV„ 


r c 

5 . 1 . 


dx 


x 


cosh {LiyJ x 2 + (x 2 lq _ 1 ) — cos ay 


f 


x , (1 — y 2 Pi) cos (2yzj) + 2y/3j sin (2yzj) 


dy 


,2\(i- p)/2 


(1 + y 2 /3])(x 2 -y 2 ) [ 

For Dirichlet and Neumann boundary conditions we obtain 

,n(i)__ 2eL/ sin og ^ f°° xP/ 2+1 J p/2 (2xzi 


O' 


= =F- 


_II_ 1 IH_ i _ \ f dx _ - IJ/ ~ 1 ' _, 

(Ax)Pl 2 + x V q z v J 2 n 9 _ 1 Jo cosh (Li^Jx 2 + UJ 2 ^) - cos ai 


(3.31) 


(3.32) 
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In figure [fl for the simplest Kaluza-Klein model with a single compact dimension of the length L 
and with the phase a (D = 4), we have plotted the total current density, L°(j l )j/e, for a massless 
scalar held in the geometry of a single plate as a function of the distance from the plate and of the 
phase d. The left/right panel correspond to Dirichlet/Neumann boundary conditions. As has been 
already noticed before, in the Dirichlet case the total current density vanishes on the plate. 



Figure 1: The total current density, L D (j l )j/e, in the topology R 3 x S' 1 for aD = 4 massless scalar 
held with Dirichlet (left panel) and Neumann (right panel) boundary conditions in the geometry of a 
single plate, as a function of the phase in the quasiperiodicity boundary condition and of the distance 
from the plate. 

For the same model, hgure [2] presents the plate-induced contribution to the current density as a 
function of the distance from the plate for various values of the coefficients in the Robin boundary 
condition (left panel) and as a function of the ratio fij/L (right panel). The numbers near the curves 
on the right panel correspond to the value of fdj/L. The left panel is plotted for the hxed value of 
the relative distance from the plate Zj/L = 0.3. On both panels, the dashed curves are plotted for 
Dirichlet and Neumann boundary conditions. For the phase in the quasiperiodicity condition we have 
taken d = 7 t / 2 . On the right panel, for the values of fij/L between the ordinate axis and the vertical 
dotted line (Pj/L = 1/d) the vacuum is unstable. 


4 Current density between two plates 


Now we turn to the geometry of two plates. In the region a\ ^ x p+l ^ < 22 , by using the formula (12.2811 
for the Hadamard function, the VEV of the current density is decomposed as 




dy 


(y 2 -^ l q ) {p 1 )/ 2 g(zj,iy) 

ci(ay)c2(ay)e 2a y — 1 


(4.1) 


Here, the second term in the right-hand side is induced by the plate at x p+1 = ay, j' / j. 

Extracting from the second term in the right-hand side of (14.11) the part induced by the second 
plate when the first one is absent, the current density is written in a more symmetric form: 

(/} = </>0 + L >'" + A(A (4.2) 

3 = 1,2 
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Figure 2: The plate-induced contribution to the current density for the model corresponding to figure 
|T]as a function of the distance from the plate (left panel) for different values of the ratio /3j /L (numbers 
near the curves) and as a function of (3j/L (right panel) for Zj/L = 0.3. The dashed curves correspond 
to Dirichlet and Neumann boundary conditions and the graphs are plotted for a = 7r/2. 


where the interference part is given by the expression 


eC, 


A (j l ) = J2 kl I d y (y 2 - W n,) 


Eni 2 + Ylj = 1,2 

O J ’ 


=-2 yzj 


/ c j( a v) 


ci(ay)c2(ay)e 2a y — 1 


(4.3) 


By taking into account the expression for the current density in the geometry of a single plate, for the 
total current density we can also write 


( f) = (j l )o + ^rr^2 k i [ dy (y 2 - ul q Y 2 

- v q nq 

2 + J2j=i,2 c j( a y) e2yZj 


ci(ay)c2(ay)e 2ay — 1 


(4.4) 


For special cases of Dirichlet and Neumann boundary conditions on both plates the general formula 
is simplified to 


(j l ) - (j l )o + 


/•OO 

^2 k i / dy(y 2 -ul q 


2 "F Ej=i ,2 e2yZj 

e 2ay _ l 


(4.5) 


where, as before, the upper and lower signs correspond to Dirichlet and Neumann boundary conditions, 
respectively. In particular, for Dirichlet boundary condition the part induced by the second plate 
vanishes on the first plate. Note that in the system of two fields with Dirichlet and Neumann conditions 
the distribution of the total current density in the region between the plates is uniform and the current 
density vanishes in the regions z < a\ and z > ci 2 - Another form for (14.511 is obtained by making use 
of the expansion 


1 


e 2 ay _ i 


£■ 

n =1 


—2 nay 


After the integration over y we get 

2e/V q 


(?) = 0 >o + 


(2it)p/ 2+ 1 ^ A 

v ' n= 1 n, 


£ ^2 ^ w S 9 [2/e (2nauj aq ) =F f$(2(na - Zj)u nq )]. 


(4.6) 


(4.7) 


.7=1,2 
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A similar representation for the interference part A (j l ) is obtained from (14.71) by the replacement 


—Zj. For Dirichlet boundary condition, on the plates, z = dj, one has 


AOVcu = 


2e/V q 

( 2 tt)p/ 2+ 1 


^2k^ q fR(2au nq ). 


(4.8) 


Combining this result with the formulas for single plates, we see that in the case of Dirichlet boundary 
condition the total current vanishes on the plates: (j l ) z=a . = 0 . 

An equivalent representation for the current density in the region between the plates and for Robin 
conditions is obtained by using the representation (1A.6I) for the corresponding Hadamard function: 

/ -i\ — / -l\ I 2 1 p ^ 2 eLi sin (nd;) f 00 , 

3 ~ {l)i y 


Wp/ 2+1 {nLi y/y 2 -!*}^) 


g(zj,iy). 


ci ( ay)c 2 (ay)e 2ay - 1 a V ~ J ’ 

Combining the expressions (13.281) and (14.91) . for the total current density we find 


(4.9) 


/ -l\ _ /-l\ i 2 p/,2 eLi \ - sin (no:;) \ - 

U / - >o t- ^p/2+iy z ^ (ULAP+ 1 ^ 

y n=l v ' n a _ 

2 + E#=i ,2 e 2yz iCj(ay) 


dy 




'—^-Wp^+iinLiy/y 2 ~ 


ci{ay)c 2 (ay)e 2a y ~1 ^ prz+LV ^ l M y ^ 4 - 10 ) 

Now, by taking into account the expression (|3.28l) for the single plate induced part, from (14.91) for the 
interference part we get 


i 2 v ! 2 eLi sin (non) v" f°° r , 

A 0) = T774XT77 E tTTTF+I E / ^ 


^p/2+ly (nL l )P+ 1 

y n =1 7 n„_ 




. 2 + E = 1 . 


=- 2 y 2 :j 


Ah (ay) 


- 1 ' w -i^ n H y2 ~ 


( 4 . 11 ) 


The equivalence of the representations (I4.4j) and (14.91) can be seen directly by using the formula (|3.30l) 
in a way similar to that for the geometry of a single plate. 

For Dirichlet and Neumann conditions, after using the expansion (14.61) . the integral over y in ()4.10p 
is expressed in terms of the MacDonald function and one gets the representation 

9 (l-p)/2 r2 °° 

^ = E+W E nsin ( n ^) E <-i 

° 1 n 9 _i 


1 n= 1 


x ( Wn 8 -i V 4 ( ra ) 2 +n2L ?) 

r=—oo ^ 

=F/p+ 3 (wn,.! 4(ra - 2 + ai) 2 + n 2 L 2 ) j 


( 4 . 12 ) 


where we have taken into account the expression (13.31) for the current density in the boundary-free 
geometry. In the model with a single compact dimension with the length L and for a massless field, 
from (14.121) we find 

i-i\ 2r((D + l)/2)e y- y-^ . , 

0) = ^(D+D/2 L D E E 


71=1 r = — C 50 


f £?+l D + l 1 

x < [4(ra/L ) 2 + n 2 ] 2 =F [4(ra — z + a\) 2 /L 2 + n 2 ] 2 >. (4.13) 
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In the case of Dirichlet boundary condition on the left plate, x p+l = aq, and Neumann boundary 
condition on the right one, x p+l = a 2 , the corresponding formulas are obtained from (14.121) and (14.131) 
with the upper sign, adding the factor (— l) r in the summation over r. The corresponding current 
density vanishes on the left plate. From (14.121) we can also see that the normal derivative of the current 
density vanishes on the plates for both Dirichlet and Neumann boundary conditions. 

In the limit a <C Li, i 7 ^ l, the dominant contribution to the series over n g _i in (14.111) comes from 
large values of \rii\, i 7 ^ l, and we can replace the summation by the integration in accordance with 

9 (' 4 7r )( 1- '?)/ 2 y r °° _ 

Y /("-,-■) -> r (( 9 - 1)/2)l! 1 *“^ 2 /(V » 2 + ™ 2 ). (4.14) 

n o —1 


Changing the integration variable y to x = \Jy 2 — u 2 , we introduce polar coordinates in the (u, x)- 
plane. After the integration over the polar angle, we get 


A{j l ) « A(j l ) R D xS i, 


(4.15) 


where A(j l ) R D xR i is the corresponding quantity in the geometry of a single compact dimension with 
the length L/. The expression for A(j l ) r d xS i is obtained from (14.111) taking p = D — 2, V q = Li, 
Lo nq _j = m, and omitting the summation over n g _ 1 . If, in addition, am <C 1, one finds 


A (j l 


2e 


V- sin (nai) f 

(nL t )D-i J 0 


(2vr 


1^/2 


(nL?) 


dy 


2 + E j= i, 2 e- 2 ^ /a /cj(y) 

cily)c 2 {y)e 2 y - 1 


w D /2 (nLiy/a). 


(4.16) 


Now let us also assume that a <C Li, m _1 , for alH = p + 2,..., D. This means that the separation 
between the plates is smaller than all other length scales in the problem. In order to estimate the 
integral in (14.161) . we note that for a fixed b and for A —>• + 00 , the dominant contribution to the 
integral / 0 °° dy f{y)e^ by w R ,/ 2 {\y) comes from the region with y < a/L. By taking into account that 


roo 

/ dye~ by w D/2 {\y) 

■Jo 


2 d / 2 X d T((D + l)/2) 
0F(6 2 + A 2 ) (d+1)/2 ’ 


to the leading order we get 


/»oo 

/ dy f(y)e~ by w D/2 (Xy) 

Jo 


2 d/2 

-r=-r((D + l)/2)/(0). 

V^A 


For the integral in (14.161) we take 6 = 2 and 


f(y) 


2 + E j= i , 2 e- 2 ^ /a / Cj -(y) 

ci(y)c 2 (y) - e~ 2y 


(4.17) 


(4.18) 


(4.19) 


In the case of non-Neumann boundary conditions one has /(0) = 1 and, hence, 


a ,J\ _ 2 er((D + l)/ 2 ) ^ sin (non) 
U ^P+i)/2 L d 


(4.20) 


Combining this result with the expressions from the previous section for the geometry of a single plate, 
we conclude that lim a _>o(y) = 0; he., for non-Neumann boundary conditions the total current density 
in the region between the plates tends to zero for small separations between the plates. For non- 
Neumann boundary condition on one plate and Neumann boundary condition on the other we have 
/(0) = —1 and the corresponding formula is obtained from (14.201) changing the sign of the right-hand 
side. In this case we have again lim a _ ) .o(j i ) = 0. 
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For Neumann boundary condition on both plates, for the function in (14.191) we have f(y) 2 /y, 
y —^ 0. In order to obtain the leading term in the asymptotic expansion for small values of a it is more 
convenient to use the expression (14.131) with the lower sign instead of the right-hand side of (14.161) . 
For small a/L the dominant contribution in ()4. 131) comes from large values of r and, to the leading 
order, we replace the corresponding summation by the integration. For the leading term this gives 


</> 


2eT(D/2) y 
7TD/2 L D-l a 2^ 


sin [na) 
n D ~ l 


(4.21) 


and for Neumann boundary condition the current density diverges in the limit a —>• 0 like 1/a. The 
described features in the behavior of the vacuum current density, L D (j l )/e , in the region between the 
plates located at z = 0 and z = a, as a function of the separation between the plates, is illustrated 
in figure [3] for a D = 4 massless scalar field in the model with a single compact dimension of the 
length L and of the phase a. The graphs are plotted for z = a/2 and a = n/2, in the cases of 
Dierichlet (D), Neumann (N) boundary conditions on both plates, for Dirichlet boundary condition 
at z = 0 and Neumann boundary condition at z = a (DN), and for Robin boundary conditions with 
fij/L = —0.5 and Pj/L = —1 (numbers near the curves). At large separations between the plates, the 
boundary-induced effects are small and the current density coincides with that in the boundary-free 
geometry. 



Figure 3: The VEV of the current density in the region between the plates evaluated at z = a/2, as 
a function of the separation between the plates. The graphs are plotted for Dirchlet and Neumann 
boundary conditions on both plates, for Dirichlet condition on the left plate and Neumann condition 
on the right one, and for Robin boundary conditions with the values of fij/L given near the curves. 
For the phase we have taken the value a = ir/2. 

In figure [U in the model with a single compact dimension of the length L and for a D = 4 
massless scalar field with Dirichlet (left panel) and Neumann (right panel) boundary conditions, we 
have plotted the total current density as a function of the ratio z/a in the region between the plates. 
The numbers near the curves correspond to the values of a/L and the graphs are plotted for a = ir/2. 
The features, obtained before on the base of asymptotic analysis, are clearly seen from the graphs: the 
current density for Dirichlet/Neumann scalar decreases/increases with decreasing separation between 
the plates and for Dirichlet scalar it vanishes on the plates. 

The same graphs for Dirichlet boundary condition on the left plate and Neumann condition on the 
right one are presented on the left panel of figure [5j The right panel in figure [5] is plotted for Robin 
boundary condition on both plates with (3\/L = fo/L = — 1. In the Robin case, the current density 
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Figure 4: The current density between the plates as a function of the relative distance from the left 
plate in the model with a single compact dimension. The graphs are plotted for a massless field 
with the parameter a = 7r/2 and with Dirichlet (left panel) and Neumann (right panel) boundary 
conditions. The numbers near the curves correspond to the values of a/L. 


decreases with the further decrease of the separation between the plates and it tends to zero in the 
limit a —>• 0, in accordance with the general analysis described above. 




Figure 5: The same as in figure [4] for Dirichlet boundary condition on the left plate and Neumann 
condition on the right one (left panel). The right panel is plotted for Robin boundary condition on 
both plates with f3\/L = fa/L = —1. 


5 Conclusion 

In the present paper we have investigated the influence of parallel flat boundaries on the VEV of 
the current density for a charged scalar held in a hat spacetime with toroidally compactihed spatial 
dimensions, assuming the presence of a constant gauge held. The effect of the latter on the current 
is similar to the Aharonov-Bohm effect and is caused by the nontrivial topology of the background 
space. Along compact dimensions we have considered quasiperiodicity conditions with general phases. 
The special cases of twisted and untwisted helds are the configurations most frequently discussed in 
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the literature. By a gauge transformation, the problem with a constant gauge held is mapped to the 
one with zero held, shifting the phases in the periodicity conditions by an amount proportional to 
the magnetic flux enclosed by a compact dimension in the initial representation of the model. On the 
plates we employed Robin boundary conditions, in general, with different coefficients on the left and 
right plates. The Robin boundary conditions for bulk helds naturally arise in braneworld scenario and 
the boundaries considered here may serve as a simple model for the branes. 

We considered a free held theory and all the information on the properties of the vacuum state 
is encoded in two-point functions. Here we chose the Hadamard function. The VEV of the current 
density is obtained from this function in the coincidence limit by using (12.6jl . For the evaluation 
of the Hadamard function we have employed a direct summation over the complete set of modes. 
In the region between the plates the eigenvalues of the momentum component perpendicular to the 
plates are quantized by the boundary conditions on the plates and are given implicitly, in terms of 
solutions of the transcendental equation (12.171) . Depending on the values of the Robin coefficients, 
this equation may have purely imaginary solutions y = Hyi- In order to have a stable vacuum with 
((p) = 0, we assume that uq > yi. Compared to the case of the bulk with trivial topology, this 
constraint in models with compact dimensions is less restrictive. The eigenvalues of the momentum 
components along compact dimensions are quantized by the periodicity conditions and are determined 
by ()2.12|) . The application of the generalized Abel-Plana formula for the summation over the roots 
of (12.171) allowed us to extract from the Hadamard function the part corresponding to the geometry 
with a single plate and to present the second-plate-induced contribution in the form which does not 
require the explicit knowledge of the eigenmodes for fc p +i (see (12.281) '). In addition, the corresponding 
integrand decays exponentially in the upper limit. A similar representation, (I2.34|) . is obtained for the 
Hadamard function in the geometry of a single plate. The second term in the right-hand side of this 
representation is the boundary-induced contribution. An alternative representation for the Hadamard 
function, (jA. 6 j) . is obtained in Appendix, by making use of the summation formula (IA.1I) . The second 
term in the right-hand side of this representation is the contribution induced by the compactification 
of the Ith dimension. 

The VEVs of the charge density and the components of the current density along uncompact di¬ 
mensions vanish. The current density along compact dimensions is a periodic function of the magnetic 
flux with the period equal to the flux quantum. The component along the Zth compact dimension is 
an odd function of the phase d/ and an even function of the remaining phases dj, i ^ l. First we 
have considered the geometry with a single plate. The VEV of the current density is decomposed into 
the boundary-free and plate-induced parts. The boundary-free contribution was investigated in m 
and we have been mainly concerned with the plate-induced part, given by (13.101) . For special cases 
of Dirichlet and Neumann boundary conditions the corresponding expression is simplified to (|3.12l) . 
The plate-induced part has opposite signs for Dirichlet and Neumann conditions. At distances from 
the plate larger than the lengths of compact dimensions the asymptotic is described by (13.131) and the 
plate-induced contribution is exponentially small. For the investigation of the near-plate asymptotic 
of the current density it is more convenient to use the representation (13.171) for the general Robin 
case and ()3.20[) for Dirichlet and Neumann conditions. From these representations it follows that the 
current density is finite on the plate. This property is in sharp contrast with the behavior of the VEVs 
of the field squared and of the energy-momentum tensor which diverge on the plate. For Dirichlet 
boundary condition the current density vanishes on the plate and for Neumann condition its value 
on the plate is two times larger than the current density in the boundary-free geometry. The normal 
derivative of the current density vanishes on the plate for both Dirichlet and Neumann conditions. 
This is not the case for general Robin condition. The behavior of the plate-induced part of the current 
density along Ith dimension, in the limit when the lengths of the other compact dimensions are much 
smaller than Li, crucially depend wether the phases d,, i ^ l, are zero or not. For df ^ 0 one has 
ojqi 7 ^ 0 and the corresponding asymptotic expression is given by (13.231) . In this case the plate-induced 
contribution is exponentially suppressed. For dj = 0, i ^ l, the leading term in the asymptotic ex- 
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pansion, multiplied by V q /Li, coincides with the corresponding current density for (p + 2)-dimensional 
space with topology R p+l x S 1 . In the limit when the length of the Zth dimension is much larger than 
the other length scales of the model, the behavior of the plate-induced contribution to the current 
density is essentially different for the cases ujqi 7 ^ 0 and ujqi = 0. In the former case the leading term 
is given by (13.261) and the current density is suppressed by the factor e~ LlUJ ° l . In the second case, for 
the leading term one has the expression (13.271) and its behavior, as a function of L/, is power law. 
In both cases and for non-Neumann boundary conditions, the leading terms in the boundary-induced 
and boundary-free parts of the current density cancel each other. 

For the current density in the region between the plates we have provided various decompositions 
((HU), (14.2j) . (|4.4|) for general Robin boundary conditions and (|4.5I) . (|4.7I) . (14.121) for special cases of 
Dirichlet and Neumann conditions). In the case of Dirichlet boundary condition the total current 
vanishes on the plates. The normal derivative vanishes on the plates for both Dirichlet and Neumann 
cases. In the limit when the separation between the plates is smaller than all the length scales in the 
problem, the behavior of the current density is essentially different for non-Neumann and Neumann 
boundary conditions. In the former case, the total current density in the region between the plates 
tends to zero. For Neumann boundary condition on both plates, for small separations the total current 
density is dominated by the interference part and it diverges inversely proportional to the separation 
(see (|4.2ip ). The results of the present paper may be applied to Kaluza-Klein-type models in the 
presence of branes (for D > 3) and to planar condensed matter systems (for D = 2), described within 
the framework of an effective field theory. In particular, in the former case, the vacuum currents 
along compact dimensions generate magnetic fields in the uncompactified subspace. The boundaries 
discussed above can serve as a simple model for the edges of planar systems. 
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A Alternative representation of the Hadamard function 

In this section we derive an alternative representation for the Hadamard function which is well suited 
for the investigation of the near-plate asymptotic of the current density. The starting point is the 
representation (12.211) . We apply to the corresponding series over ni the summation formula [I4j 31] 

r\_ 00 COO 

y ^ g(k)f{\k\)= du[g(u) + g(—u)\f(u) 

Ll n; =—00 J 0 

+i f [/(*«)-/(-*«)] ^ p ( aj ) 

1/0 A=±l 

where ki is given by (12.121) . The part in the Hadamard function coming from the first term in the 
right-hand side of (IA.1D coincides with the Hadamard function for the geometry of two plates in 
D-dimensional space with topology R p+2 x T 9 ” 1 and with the lengths of the compact dimensions 
(L p+ 2 ,..., 1 , Li + 1 ,..., Ljj) (the Zth dimension is uncompactified). We will denote this function by 
G R p+ 2 xT q-i(x,x'). As a result, under the assumption f5j ^ 0, the Hadamard function is decomposed 
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as 


G[x,x') — G/jp +2 xTi-l (x, x') + 


Li 


dk„ 


iraV q J (2vr)P 


E 


nq-1 


E 


n =1 


An “I - cos [An + 27 j(A n )] sin X n 


/ 

CJ 


00 cosh(At-\/ u 2 — a;^ 2 ) e -\uAx l 

du v x 


(0 


ti — CJ 


(02 


E 

A=±l 


e uLi+i\ai _ ^ ' 


(A.2) 


where x^_ x = (5 


,p +2 


l -1 ,-vJ+l 


...,x ,x 


,... x D ), k g _i = (k p+ 2 , ..., ki-i,k + i,..., k D ), and _ fcf. 

Here, the second term in the right-hand side vanishes in the limit L; —>• 00 and is induced by the 
compactification of the Zth dimension from R 1 to 5 1 with the length L/. 

By making use of the relation 


E 

A=±l 


AuArd 


e uLi+i\ai _ ]_ 


2 u ^2 h r (u, Ax 1 ), 


r—1 


with 


g—ruLi , >. 

h r (Ax l ,u) =-cosh ( uAx 1 + iron ) 

u V / 


(A.3) 


(A.4) 


we rewrite the formula (1A.2I) in the form 


2Li 


g(x,x') — G RP +2 XT g-i(x,x') + _ 2 


dk r . 


^ J (2vr )p 

/*oo 

V / dy cosh(yAf)e^' Ax ^ k s-i Ax *-i 
. , */0 


nq -1 


“ A„g(z,z', X n /a)h r (Ax l , yjXI/a 2 + y 2 + 

/ 


n=l 


An + cos [An H - ^7? (An)] sin A^ 


(A-5) 


with w Pin(3 _ 1 = y/k^ + w^_ 1 . Now, by using the summation formula (12.2511 for the series over n we 
get the final representation 


G{x,x') = 


e ikp-Ax p +ik^_ 1 .Axj ! _ 1 


Gzjp+i x r*-i ( x ’ x ') + ^flf X2 1 

q r =1 J v y n,_i 

xj dy cosh(Aty)| J dugj(z, z', u)h r {Ax 1 , yju 2 + y 2 + ng _ 1 ) 

+/; 


gj(z,z',iu ) 


'2v 2+ < 


ci (cm) C 2 (ait) 


j e 2an _ ! E *^r(As , ,y« 2 -g 2 -w2 nfl _ i )j.( A . 6 ) 


In this expression, the part with the first term in the figure braces is the contribution to the Hadamard 
function induced by the compactification of the Zth dimension for the geometry of a single plate at 
x p+1 = aj and the part with the second term in the figure braces is induced by the second plate. Note 
that the contribution of the first term in the right-hand side of (jA.611 to current density along the Ith. 
dimension vanishes. 

In deriving the representation (1A.6I) we have assumed that f3j ^ 0. For this case, in the region 
between the plates, all the eigenvalues for the momentum k p+ 1 are real and in the geometry of a single 
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plate there are no bound states. For /3j > 0, in the application of the summation formula (|2.25l) to 
the series over n in ()A.5I) the contribution from the poles ±i/bj should be added to the right-hand 
side of (12.251) . This contribution comes from the bound state in the geometry of a single plate at 

x p+1 = ctj. For this bound state the mode function has the form (x) ~ e~ Zj ^ j e lk ll' x ll =F * a; k f with 
ui^' = kp T — 1 //3j . Assuming that ujqi > 1 //3j, the contribution from the bound state to the 
Hadamard function in the geometry of a single plate is given by the expression 

qM( t r '\ - -\z+z'-2a i \/0 i Y' [ dk P V [°° H ik p -Ax p +il4_ 1 -Ax^_ 1 

bj[ ’ } ~ vVJj £[J (2vr) P f ,/ 0 

r ~ i n ,-i 

x cosh(xAt)h r (Ax l , \Jx 2 + k£ + — 1//3?), (A.7) 

where 6{x) is the Heaviside unit step function. In the case ojqi < 1 //3j < uo the corresponding 
expression is more complicated. 
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